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Using a constrained formalism for Einstein equations in Dirac gauge, we propose to 
' compute excised quasistationary initial data for black hole spacetimes in full general 

relativity. Vacuum spacetime settings are numerically constructed by using the isolated 
horizon formalism; we especially tackle the conformal metric part of our equations, as- 
suming global stationarity. We show that a no-boundary treatment can be used on the 
horizon for the equation related to the conformal metric. We relate this finding to previ- 
ous suggestions in the literature, and use our results to assess the widely used conformally 
flat approximation for computing black hole initial data. 

> 

o . 

CN| • 1. Describing stationary black hole data with isolated horizons 

o ■ 

■ Aiming at the computation of black hole initial data in numerical relativity, this 

work suggests the use of an excision approach, by solving Einstein equations outside 
two-spheres assumed to encircle black hole singularities. In the spirit of the mem- 
brane paradigm^ we describe black holes as physical objects characterized by their 
horizons. We shall use here the local characterization for black hole isolated hori- 
zons, associated to a nonevolving black hole region. Isolated horizons are defined 2 
as three-dimensional null tubes foliated by marginally trapped surfaces (geomet- 
rical expansion along the outer future null normal vanishes), with the additional 
restriction that the extrinsic geometry of the tube is not evolving along the null 
generators. Those features will translate in our computation as boundary condi- 
tions on our chosen (3+1) variables associated to the metricP^ 

We choose to write Einstein equations following a (3+1) fully constrained formal- 
ism,^ with maximal slicing and spatial Dirac gauge. In the vacuum, asymptotically 
flat, single black hole spacetime setting, we end up with a subset of three elliptic 
equations (two scalar, one vector) to solve, associated to the classical lapse TV, shift 
f3 l and conformal factor tp (3+1) variables, and expressing the Einstein constraints 
as well as the gauge choice. Those equations require boundary conditions at the 
inner excised surface, provided by the prescription of this surface to correspond to 
a slice of an isolated horizon. Two free parameters allow to fix the mass and angu- 
lar momentum of the black hole. Solving for those three PDE on a spacelike slice, 
alongside with setting the conformal 3-metric ■y 1 -' = ip 4 ^ (ip being the conformal 
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factor) to be the usual flat metric leads to vacuum black hole initial data with the 
classical conformal flatness hypothesisP 

2. Solving for the conformal metric: the no-boundary approach 

Using a global stationarity hypothesis, our formalism provides a fourth elliptic-like 
equation for the deviation tensor h 1 ^ , being the non-flat part of the conformal metric 
(h^ = 7 lJ — , with /*■?' the reference flat metric) P 

Ahij ~ j^Wrf* = S ij (h^,N,^,p). 

with C the classical Lie derivative. Going beyond the conformal flatness assumption 
requires to solve this equation, and particularly to treat the boundary condition 
problem. Concerning boundary conditions for the conformal metric in excised initial 
data, a few proposals have been formulated in the literature. Jaramillo^ proposed 
to express a geometrical property of the isolated horizon, namely the invariance of 
the extrinsic geometry. Cook and Baumgarte^ expressed that this boundary value 
could be arbitrarily prescribed with no effect on the physical content of the data. 
On the contrary, we describe in this work a no-boundary treatment, arguing that 
the operator is weakly singular at the horizon, and thus can be inverted without 
additional information at the boundary. The data obtained are coherent with the 
prescription of Jaramillo, but not with the one of Cook and Baumgarte. Details on 
the resolution of this tensor equation, including an extraction of the two relevant 
scalar degrees of freedom in Dirac gaugeP can be retrieved from the article related 
to this reportP Let us note that this result has been obtained one more time since 
this studyP^l with a different approach in the treatment of the conformal metric. 

3. Results and pertinence of the study 

The data obtained from our computations verify the Einstein constraints, as well 
as the (3+1) evolution equation in stationarity, so that we are able to reconstruct 
the whole spacetime by simple time translation. The initial data settings (single 
black hole in vacuum spacetime) suggest that we should obtain from this study 
the Kerr solution in Dirac gauge (which has no analytic formulation). We test this 
assumption with several markers: first, although computations are performed with 
no spatial symmetry hypothesis, the final data arc indeed globally axisymmetric. 
Computation of the virial integral at infinity^ shows that contrary to the corre- 
sponding conformally flat data, our spacetime appears to be stationary. Finally, 
performing a source multipole decomposition 1 ^ of the isolated horizon slice, and 
comparing with values for an horizon of Kerr geometry constructed analytically, we 
show that our horizon geometry is the one of a Kerr black hole to a very reasonable 
accuracy (see Fig. [1]). Equating the isolated horizon multipoles suffices to identify 
both spacetimesP^ The corresponding data with a conformal flatness hypothesis 
show however a clear discrepancy. 
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Fig. 1. Computations of the rescaled second-order mass multipole (M2/M 3 ) (M being the ADM 
mass computed at infinity) and the rescaled third-order angular multipole (J3/M 4 ). Relative differ- 
ences with respect to values for a Kerr-Schild analytical horizon are displayed in the conformally 
flat (CF) or the non-conformally flat (NCF) setting, against the usual (dimensionless) rotation 
parameter for a Kerr spacetime. 

To our knowledge it is the first time the full 3-metric (including the conformal 
part) for black hole data is obtained numerically using only a stationarity assump- 
tion. Further numerical work using this approach to binary black hole settings could 
lead to improvements in the computation of initial data for such systems. 
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